Mon. Not. R. Astron. Soc. 000, ITHTOl (20071 Printed 15 January 2013 (MN WT^i style file v2.2) 



o 

(N 



Gamma-rays from nebulae around binary systems 
containing energetic rotation powered pulsars 

W. Bednarek 1 & J. Sitarek 2 

1 Department of Astrophysics, The University of Lodz, ul. Pomorska 149/153, 90-236 Lodz, Poland, bednar@astro.phys.uni.lodz.pl 
2 IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain, jsitarek@ifae.es 



Accepted . Received ; in original form 



ABSTRACT 

We consider nebulae which are created around binary systems containing rotation 
powered pulsars and companion stars with strong stellar winds. It is proposed that 
the stellar and pulsar winds have to mix at some distance from the binary system, 
defined by the orbital period of the companion stars and the velocity of the stellar 
wind. The mixed pulsar-stellar wind expands with a specific velocity determined by 
the pulsar power and the mass loss rate of the companion star. Relativistic particles, 
either from the inner pulsar magnetosphere and/or accelerated at the shocks between 
stellar and pulsar winds, are expected to be captured and isotropized in the reference 
frame of the mixed wind. Therefore, they can efficiently comptonize stellar radiation 
producing GeV-TeV 7-rays in the inverse Compton process. We calculate the 7-ray 
spectra expected in such scenario for the two example binary systems: J1816+4510 
which is the redback type millisecond binary and LS 5039 which is supposed to contain 
energetic pulsar. It is concluded that the steady TeV 7-ray emission from J1816+4510 
should be on the 100 hr sensitivity limit of the planned Cherenkov Telescope Array, 
provided that e ~ 10% of the rotational energy lost by the pulsar is transferred to 
TeV electrons. On the other hand, the comparison of the predicted steady TeV 7- 
ray emission, expected from 7-ray binary LS 5039, with the observations of the TeV 
emission in a low state, reported by the H.E.S.S. Collaboration, allows us to put 
stringent upper limit on the product of the part of the hemisphere in which the mixed 
pulsar-stellar wind is confined, A m ; x , and the energy conversion efficiency, e, from the 
supposed pulsar to the TeV electrons injected in this system, A mix -e < 1%. This lower 
limit can be understood provided that either the acceleration efficiency of electrons 
is rather low (e ~ 1%), or the parameters of the stellar wind from the companion 
star are less extreme than expected, or the injection/ acceleration process of electrons 
occurs highly anisotropically with the predominance towards the companion star. 

Key words: binaries: general — pulsars: general — radiation mechanisms: non- 
thermal — gamma-rays: theory 



1 INTRODUCTION 

Some binary systems contain rotation powered pulsars 
which are energetic enough to prevent accretion of matter 
from stellar winds of companion stars. These pulsars can 
create strong pulsar winds which collide with the winds of 
the companion stars. As a result, a shock structure appears 
within the binary system on which particles can be accel- 
erated to relativistic energies. A few types of such binaries 
have been discovered up to now. The most famous is a class 
of " Black Widow" binaries in which energy output from the 
millisecond pulsar can even partially evaporate the compan- 
ion star. The companion stars in these systems are char- 
acterised by small masses, below ~ O.IM©, and relatively 



strong winds induced by energy released by the pulsar. In 
another type of such binaries, the so-called "Redback" bi- 
nary systems, the energetic millisecond pulsa r forms a binar y 
system with a low mass ~ 0.2 M star (see lRobertsll201lf ). 
Up to now, only one binary system has been discovered in 
the Galaxy to contain a classical 7-ray pulsar on a close or- 
bit ar ound the massive, Be t ype, s tar (i.e. PSR 1259-63/SS 
2883, Ijohnston et~ai1 1 19921 . \19d4) . It has been proposed 
that other binary systems c an also contain ejecting cl assical 
pulsa rs, e.g. LSI 303 +61 dMaraschi fc Trevesl ll98l|V Cy g 



X-3 (|Harding fc Gaisserl Il990l ) or LS 5039 (|Dubusl I2OO6I ) . 

These binaries emit 7-rays in the GeV-TeV energy range, 
clearly modulated with the period of t he binary system 
l| Albert et al. 2009; Ahar onian et alfeOOrJ ) . In the most pop- 
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ular scenario, these 7-rays are produced by relativistic elec- 
trons which up-sc atter stellar radiation in the Inverse Comp- 
ton process, (e.g. iMaraschi fc Treveslll98ll ; pTavani fc Aronsl 
1 19971 ; lBednarek| [l997l ). Since the stellar radiation field is 
anisotropic for electrons ejected from the pulsar or acceler- 
ated at the shock, the emission has to be modulated with the 
period of the binary system. However, a part of relativistic 
electrons in the pulsar wind can escape from the binary sys- 
tem without significant interactions, i.e. those moving close 
to the outward direction. In the case of low mass millisecond 
pulsar binaries, most of the accelerated electrons can escape 
radially from the vicinity of the companion star, due to rela- 
tively weak radiation field. Farther from the binary systems, 
the winds from the pulsar and the companion star have to 
mix due to the rotation of the binary system. As a result, 
it is expected that such mixed pulsar-stellar winds move 
together with some specific velocity determined by the en- 
ergy output from the pulsar and the material content of the 
stellar wind. Therefore, nebulae around binary systems con- 
taining rotation powered pulsars should have different struc- 
tures than nebulae around isolated pulsars. We propose that 
such binary systems should be surrounded by nebulae com- 
posed of the mixture of the pulsar and stellar winds which 
expands with non-relativistic velocity. Relativistic particles 
are expected to be captured and isotropized in the reference 
frame of this wind. 

In this paper we consider radiation processes occurring 
in such nebulae. We calculate the 7-ray spectra produced by 
relativistic electrons in the Inverse Compton (IC) scattering 
of thermal radiation from the companion star. In contrast to 
the modulated 7-ray emission from the interior of the binary 
system, we predict the presence of a steady component of 
the 7-ray emission produced in the nebula around the binary 
system. Discovery of such additional component by the fu- 
ture Cherenkov telesco pes (e.g. Cherenkov Telescope Array, 
CTA. lActis et al1l201lT ) from the binaries containing pulsars 
should allow to derive important parameters determining 
the acceleration processes in such binaries and properties of 
the emission region. The comparison of the predicted high 
energy emission with the available observations of the TeV 
7-ray emission in the low state from already known 7-ray 
binaries should allow to provide additional constraints on 
the acceleration efficiency of particles occurring around the 
rotation powered pulsars in the binary systems. 



2 NEBULAE AROUND BINARY SYSTEMS 

As we have outlined above we consider specific type of nebu- 
lae around binary systems containing rotation powered pul- 
sars and companion stars with substantial mass loss rate. 
In fact, in such type of binaries a double shock structure 
appears already within the binary system in collision of the 
pulsar and stellar winds. Its parameters are determined by 
the energy loss rate of the pulsar and the power of the 
stellar wind. In the case of spherical winds, the distance 
of the shock from the star is defined by the parameter, 
v = ( L m i/c)/M v v v « 0.05L 3 6/(M-7«8) (|Girard fc Willsonl 
1 19871 ). where the pulsar wind power is L pu i = 10 Lza 
erg/s, the stellar wind mass loss rate is M w = 10 -7 M_7 
M©/yr and its velocity is u w = 10 8 i>8 cm/s. The closest, 
radial distance of the shock from the pulsar is then given 



by p — Dby/rj/ (1 + -Jrf), where Db is the separation of the 
binary stars. For the example scaling parameters applied 
above, the shock should be bound around the pulsar. In 
the case of binary systems containing energetic millisecond 
pulsars, the shocks are usually bound around the low mass 
companion stars due to relatively small power of the stel- 
lar winds in respect to power of the winds around massive 
stars. These shock structures can be much more complicated 
in the case of no n-spherical stellar winds of massive star s of 
the Be type fsee lSierpowska-Bartosik fc Bednarekllioosl ). In 
this last case the location of the shock structure can change 
drastically depending on the location of the pulsar within 
the polar or equatorial stellar wind. 

The shock structure tends asymptotically to a cone 
characterised by a half-opening angle ip — 2.1(1 — 
. T 7m 5 /4)r?m 3 rad, where r? m = min(r;,7) _1 ) (see lBall fc Doddl 
l200ll ; lEichler fc Usovlll993h . However, due to the rotation 
of the binary system, the cone should be bound at large 
distances from the stars (e.g. Dubus 2006, Bosch-Ramon et 
al. 2012) provided that the velocity of the stellar wind is 
large in respect to the rotational velocity of the companion 
star. For r\ > 1 the pulsar wind completely overtakes the 
stellar wind, creating a torus like structure in the equatorial 
plane of the binary system (see. Fig. 1 for schematic geome- 
try). In this case, the pulsar wind can only propagate freely 
along the polar regions of the binary system. In the opposite 
case r] < 1, the stellar wind overtakes the pulsar wind. How- 
ever, again in the equatorial region the pulsar wind mixes 
with the stellar wind as in the previous case. There are sig- 
nificant differences between the mixed pulsar-stellar wind in 
the equatorial region and the freely expanding pulsar wind in 
the polar region. Relativistic electrons are isotropized in the 
reference frame of a relatively slowly expanding mixed wind. 
These electrons can interact efficiently with the thermal ra- 
diation from the nearby companion star. Relativistic elec- 
trons in the pulsar wind, expanding in polar directions, move 
rectilinear from the binary system. These electrons can be 
isotropized only far away from the binary system where the 
density of stellar photons is already very low (as in the case 
of P WNe around isolated classical pulsars) . These electrons 
are not expected to produce efficiently 7-rays by scattering 
stellar radiation. They finally scatter the Microwave Back- 
ground Radiation and the infrared radiation in the Galaxy 
producing extended 7-ray sources which are more difficult 
to detect with the Cherenkov telescopes. If the parameter 
77 is lower than unity, then the equatorial regions of the bi- 
nary system are dominated by the mixed pulsar-stellar wind 
but the polar regions are dominated by the stellar wind. In 
this case all electrons in the pulsar wind are expected to be 
isotropized at some distance from the binary system. 

We calculate the factor which determines a part of the 
hemisphere (the solid angle divided by 4tv) in which the 
mixed pulsar-stellar wind (the toroidal structure in the plane 
of the binary system, see Fig. 1) is confined. This factor is 
calculated for the known value of the half opening angle, ip, 
of the shock created within the binary system in collisions 
of these two winds. This part of the hemisphere is equal to 
A m i x = Af2 m i x /47r = (1 — cost/)), where Af2 m i x is the solid 
angle in which the mixed pulsar-stellar wind is confined. 
Both winds are expected to mix completely at the distance 
from the binary system which depend on the velocity of the 
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Figure 1. Schematic representation (not to scale) of the vicinity of the binary system containing an energetic pulsar and a companion 
star (seen from the reference frame of the pulsar). The top view (i.e. from the pole of the binary system) is shown on the left and the 
back view (from the side of the binary plane) is shown on the right. On the left: The dotted line shows the orbit of the star. The pulsar 
creates energetic wind which collides with the stellar wind. As a result, a double shock structure is created (marked by the dashed curve). 
The stellar (pulsar) wind creates a cometary tail which interacts with the pulsar (stellar) wind. At some distance from the binary system 
(signed by D m i x ) both winds mix together. The mixed stellar-pulsar wind propagates with the common velocity estimated from Eq. [3] 
Rclativistic particles accelerated by the pulsar move rectilinear in the pulsar wind up to the shock structure. In this region many of 
them are not able to interact with the stellar radiation due to the small interaction angles. However, in the mixed stellar-pulsar wind 
they are isotropized. They start to lose efficiently energy on comptonization of the stellar radiation. Moreover some particles can be also 
accelerated at the shock between the pulsar and stellar winds in the Fermi acceleration process. On the right: This same picture but 
viewed from the plane of the binary system for the case of the parameter r\ > 1 . In this case the pulsar wind can escape freely only along 
the polar regions of the binary system. The stellar and pulsar winds have to mix at the distance D m i x in the region of the plane of the 
binary system. 



stellar wind and the rotational period of the binary system, 

-Dmix ~ T or b«w ~ 8.6 x 10 12 r d w g cm, (1) 

where the period of the binary system is r or b = lrd days. In 
the above formula we assumed that the stellar wind velocity 
is much greater than orbital velocity of the star which is 
usually the case for the considered binaries. As a result of 
the mixing process, the pulsar wind becomes uploaded with 
the stellar wind matter. The mixture of both winds starts 
to move together with the common velocity, w; onv , estimated 
from the energy conservation, 

Lpul + isw = MWcnv/2, (2) 

where L sw = M w u^,/2 is the kinetic power of the stellar 
wind. Here we neglect the energy loss of the bulk motion of 
the mixed winds on the isotropization of the direction of the 
particles in it. If the stellar wind power dominates over the 
pulsar wind power, then iaj c1w ~ v w . In the opposite case, 

w env = ^2L pnl /M « 5.4 x 10 8 (L 36 /M-7) 1/2 cm/s, (3) 

resulting in a non-relativistic bulk motion of the mixed 
winds for the assumed above example parameters. 

Relativistic electrons accelerated in the inner binary 
system, either by the pulsar itself or on the shock formed 
in the pulsar and stellar winds collision, are immersed in 
the mixed stellar-pulsar wind. We estimate the maximum 
energies of the electrons for which they are captured and 
confined by the mixed wind magnetic field. It is assumed 



that the magnetic field in the mixed wind is determined by 
the magnetic field in the stellar wind. It has a dipole struc- 
ture only very close to the stellar surface. Farther out, it has 
a radial structure and finally reaches a toroi dal structure due 
to th e rotation of the star (for details see lUsov fc Melrose! 
Il992i ). The magnetic field strength with the toroidal struc- 
ture can be approximated as a function of the distance from 
the star by the following formula, 

B(D) « O-OlB2.Ru/-D13 G, (4) 

assuming that the transition between the region of the 
toroidal and the radial magnetic field occurs at the dis- 
tance of ~ 10 stellar radii (where the radius of the star 
is R* = lO^-Rn cm). In this formula, the surface magnetic 
field of the star is B, = IOOB2 G, and the distance from the 
star is expressed as D = 10 13 Bi3 cm. 

We assume that relativistic electrons are captured by 
the wind (and advected with the wind) if their Larmor 
radius, Rl, is lower than the characteristic distance scale 
given by the distance from the massive star, i.e. _Rl ~ 
3 x 10 9 B T cv/B(D> mix )cm < D mbc , where E = lB Tc v TeV is 
the energy of an electron, and B(D' m i x ) (measured in Gauss) 
is the local magnetic field strength in the mixed wind. Based 
on this condition, we conclude that electrons with energies, 

E < 30B 2 TeV, (5) 

should be captured in the mixed stellar-pulsar wind. 

Electrons, captured in the magnetic field of the mixed 
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stellar-pulsar wind, are isotropized in the reference frame 
of the relatively slow wind. They lose energy on the syn- 
chrotron radiation and on the IC scattering of the stellar 
radiation. The importance of the synchrotron energy losses 
can be evaluated by comparing the synchrotron cooling time 
scale with the advection time scale with the velocity of the 
winds. The advection time scale can be estimated from, 

Tadv = D/w cnv = 10 5 Di 3 /w 8 s, (6) 

where v cnv — 10 8 u>8 cm/s, and the synchrotron time scale 
from, 

r syn = E c /E syn « 290/(B 2 £) s, (7) 

where E syn = (4/3)cU B a T E 2 /m 2 w 3.5 x W^ 3 B 2 E 2 TeV/s, 
m e — 511 keV and electron energy E is expressed in TeV. 
By comparing Eq. [6] and [7] we get the limit on the electron 
energy above which the synchrotron energy loss time scale 
is shorter than the advection time scale, 

EZl « 2<dw a D l3 /{B 2 Rii) 2 TeV. (8) 

This critical electron energy is equal to E c « 
130r d vs(L36/M- 7 ) 1/2 /(B 2 Rii) 2 TeV, for the mixture 
distance of the winds given by Eq. [1] and the velocity of 
the mixed wind given by Eq. We conclude that electrons 
with TeV energies are not able to lose efficiently energy on 
the synchrotron process during their advection with the 
mixed pulsar-stellar wind. However, at larger energies the 
synchrotron energy losses should be taken into account. 

Relativistic electrons, after isotropization by the mag- 
netic field in the mixed pulsar-stellar wind, should lose effi- 
ciently energy also on the IC scattering of the stellar radi- 
ation. We estimate the optical depth for the IC scattering 
process in the Thomson regime for electrons at the distance 
of Dmix on, 

Tic/t = n*o- T c(D mi x/w nv) ~ 4.7T 4 3 7ii 1 /(rdt;8ii'8), (9) 

where the density of stellar photons is n* = 2.75 x 
10 9 T 4 3 (i?ii/(r d V8) 2 cm , and -D m ix/w e nv is the characteris- 
tic time scale spend by electrons at the distance D m i x from 
the binary system. Note that the scattering in the Thom- 
son regime is valid only for electrons with energies below 
m 2 e 4 /3fc B T ~ 100/T 4 GeV. However, in the case of LS 5039, 
the optical depth for the ~TeV electrons should be still close 
to unity, keeping in mind that the Klein-Nishina cross sec- 
tion drops inversely proportional with the electron energy. 
We conclude that electrons with TeV energies should ef- 
ficiently lose energy on the IC scattering of stellar radia- 
tion after being isotropized in the pulsar-stellar wind nebula. 
However, in the massive binary systems of the LS 5039 type, 
electrons, moving rectilinear in the pulsar wind in the out- 
ward direction from the companion star, are not expected 
to lose significantly energy already in the binary system due 
to the small scattering angles between electrons and stellar 
photons. The 7-rays produced at the distance D m i x from 
the star should escape from stellar radiation field without 
significan t absorption (e.g . see the scaling formula given by 
Eq. 26 in lBednarekl (|2009T ) for the optical depths in the stel- 
lar radiation field). 

The nebulae formed by the mixed pulsar-stellar winds 
can have very complicated structures in the case of binary 
systems in which the wind of the companion star is non- 
spherically-symmetric or the orbit has a large eccentricity. 



For example, in the case of the eccentricity not far from unity 
the distance, D m ix, differs significantly for the periastron 
and apastron passages. These distances can be estimated 
from, -DJ^ ~ TpetVw, for the periastron passage and from 
^mix ~ Tapo^w for the apastron, where r pcr and T apo are 
the time scales describing the passage of the compact object 
through a half of the orbit centered on the periastron and 
apastron, respectively. In such eccentric binaries, nebulae 
should have the shape corresponding to the shape of the 
orbit of the companion stars. 

In the present paper we outline only qualitatively the 
scenario describing the process and properties of the mix- 
ing of stellar and pulsar winds and concentrate on the ex- 
pected high energy radiation produced in such scenario. The 
details of the mixing process of the winds are very compli- 
cated. However, the first results of hydrodynamic numerical 
simulations of collisions of stellar winds are consistent with 
general features of the scenario described in this section (see 
e.g. Pittard 2009, Lamberts et al. 2012). 



3 NEBULAE AROUND SPECIFIC BINARY 
SYSTEMS 

As an example, we estimate the parameters of nebulae 
around a few binary systems containing (or expected to 
contain) energetic pulsars, i.e. classical radio pulsar PSR 
1259-63, millisecond pulsar B1957+20 in the Black Widow 
type binary system, millisecond pulsar J1816+4510 in the 
Redback type binary system, and supposed classical pulsar 
in the binary system LS 5039 belonging to the class of the 
TeV 7-ray binaries. The basic parameters of these binaries, 
collected or proposed in different works, are summarised in 
Table 1. 

Considered above binary systems are expected to be 
surrounded by nebulae with quite different parameters. We 
estimate the basic parameters of these nebulae based on the 
formulae derived in Sect. 2 (see Table 2). The best con- 
ditions for 7-ray production seems to be provided by the 
binary system of the redback type, containing millisecond 
pulsar J1816+4510, and the supposed pulsar within the bi- 
nary system LS 5039 from which already GeV and TeV 7- 
ray emission has been detected. The observed emission from 
LS 5039 is modulated with the binary period. Therefore, 
it is expected to be originated within the binary system. 
We propose that additional steady 7-ray emission compo- 
nent should be also produced in the nebula surrounding this 
system provided that it contains energetic pulsar. We cal- 
culate the steady 7-ray fluxes from the nebulae surrounding 
two binary systems mentioned above selecting them as the 
best candidates for production of observable steady 7-ray 
emission in the surrounding nebulae, i.e. those containing 
millisecond pulsar J1816+4510 and LS 5039. 

3.1 J1816+4510 

The millisecond pulsar J1816+4510, with the period of 3.2 
ms, was recently discovered in th e survey as a part of the 
Green Bank North Celestial Cap |Stovall et alj|2012h . This 
pulsar is in the binary system with the low mass compan- 
ion (mass below O.I6M0) having an orbital period of 8.7 
hrs. The distance to this binary has been estimated on 2.4 
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Table 1. Basic parameters of stars and pulsars in binary systems: the radius of the companion star i?*, its surface temperature T*, the 
mass loss rate, M, and the velocity, u w , of the stellar wind, the magnetic field strengths at the surface of the star B+, the binary orbital 
period T or t, and the spin down luminosity of the observed (or supposed) pulsar L pu i. 



Name 


R* (cm) 


T, (K) 


M (Af /yr) 


v„ (km/s) 




(G) (?) 


T"orb (days) 


£pul (erg/s) 


B1259-63 


4.2 x 10 11 


2.7 x 10 4 


2 x 10~ 7 


2 x 10 3 


10 2 


- 10 3 


1236.7 


8.3 x 10 35 


B1957+20 


10 11 


8 x 10 3 


3 x 10- 10 


7 x 10 2 


1 - 


10 3 


0.38 


7.5 x 10 34 


J1816+4510 


7 x 10 10 


2 x 10 4 


3 x 10~ 10 


7 x 10 2 


1 - 


10 3 


0.36 


5 x 10 34 


LS 5039 


6.5 x 10 11 


3.9 x 10 4 


2.6 x 10" 7 


2.4 x 10 3 


10 2 


- 10 3 


3.9 


10 37 



kpc, based on the dispersion measure (jKaplan et al.ll2012T ). 
The pulsed 7-ray emission from J1816+4510 has been also 
detected by Fermi, implying an energy conversion efficiency 
of ~ 25% from the pulsar's spin down luminosity into j- 
rays. The pulsations have two peaks separated by ~ 0.5 in 
phase. The 7-ray spectrum is measured between 0.5 GeV 
and 5 GeV. The spectrum is flat. It can be described by a 
pure power law with a photon index —2.20 ±0.07 or a power 
law with an inde x of —2.0 ± 0.1 and an exponential cut-off 
at 7.5 ± 4.0 GeV (|Kaplan et alj|2012h . The companion star 
has a radius of O.IRq a nd the surface temp erature in the 
range 1.8 - 2.1 x 10 4 K (|Kaplan et alj|2012h . We have se- 
lected this binary as one of the interesting examples for our 
farther analysis due to the large spin down luminosity of the 
pulsar, ~ 5 x 10 34 erg s" 1 and a very hot companion, which 
is unusual for the Redback type millisecond binary system. 

The likely parameters of this binary system allows us 
to estimate the shape of the shock in this binary defined by 
the parameter r\. For r\ ~ 1.2 the shock almost divide the 
volume of the binary in two equal parts. For such geome- 
try, the mixed wind dominates almost at the whole sphere 
around the binary system, i.e. A m j x ~ 1. The energies of 
electrons accelerated in this system can reach TeV energies 
(see Table 2) for the intermediate magnetic field strengths 
from the range mentioned in Table 1. The acceleration pro- 
cess of these electrons will be considered in a more detail in 
the next section. Also the optical depths for electrons at the 
distance of D m ix are in this case of the order of a few see 
Table 2). 

3.2 LS 5039 

LS 5039 is a well known binary system containing O type 
massive sta r and a compact obj ect on the orbit with a period 
of 3.9 days dCasares et al. 20051). T he distance to the system 
is 2. 9 ±0.8 kpc (|Casares et alj|20l3 ). The vicinity of LS 5039 
shows radio morphology whic h indicates a presenc e of an en- 
ergetic pulsar in this system (|Moldon et al. 2012j) as earlier 
suggested by iDubusl (|2006t ) . This binary system has been 
detected in the TeV 7-rays jAharonian et al]|2005l ) and also 



at GeV energies by Fermi (|Abdo et al. 20091 ). showing per 



sistent emission modulated with the orbital period of the bi- 
nary system. The GeV and TeV 7-ray light curves are clearly 
anticorrelated. The 7-ray spectrum extends up to ~ 20 TeV 
indicating the presence of multi-TeV particles within this bi- 
nary system. The spectral slope changes dramatically with 
the o rbital phase between -1.85 and -2.53 (|Aharonian et all 
I2006J ). We have selected this binary due to relatively small 
eccentricity (e ~ 0.35) and the well measured modulated 
TeV 7-ray emission which should allow to put constraints 



on the acceleration rate of relativistic electrons in terms of 
our model. 

For this binary system, the shock structure bends 
around the pulsar, 77 ~ 0.08 (see Table 2). This means that 
the whole pulsar wind is overtaken by the stellar wind. The 
mixed pulsar-stellar wind propagate only in a part of the 
sphere. Electrons are expected to be accelerated to several 
TeV for the likely surface magnetic field strength of the com- 
panion star. The optical depth for electrons in the radiation 
field of the companion star (in the Thomson regime) is very 
large (of the order of ~100) at the distance of D m ix- Thus, 
efficient production of 7-rays is expected even from electrons 
with energies of the order of several TeV, which already scat- 
ter stellar radiation in the Klein-Nishina regime. 



4 RELATIVISTIC PARTICLES IN NEBULAE 

We assume that relativistic particles in the nebula surround- 
ing binary system are directly injected from the pulsar in- 
ner magnetosphere (close to mono-energetic case) or they 
are accelerated in the shocks created in pulsar-stellar wind 
collisions. The Lorentz factors of the mono-energetic parti- 
cles, 7e u , are usually considered to be of the order of a few 
times lO^ 4-6 ' . It is expected that a significant part of these 
particles, those injected in the outward direction from the 
companion star, escape to the surrounding nebula and reach 
the distance of D m i x . We estimate the maximum energies of 
particles which can be accelerated at the turbulent region at 
the distance of D m ix by comparing their acceleration time 
scale and their advection time scale from the shock region 
given by Eq. 6. The acceleration time scale is given by, 



7"acc = £c/Pacc « 100£ T cv/(X-3B) S, 



(10) 



where P acc = \cE/Rl and Rl is the Larmor radius of parti- 
cles. The acceleration process is parametrised in this case by 
the acceleration efficiency \ = 10 _3 X-3 which is assumed to 
be of the order of ~ (w env /c) 2 . Applying Eq. 4, we estimate 
the maximum energies of electrons due to their advection 
from the acceleration region on, 



EZ7 w WX-3B2R11/WS TeV. 



(11) 



The maximum energies of particles, accelerated at the shock 
at the distance of D m i x , can be additionally constrained by 
their energy losses. Therefore, we compare the acceleration 
time scale of particles with their energy loss time scale. For 
the radiation energy loss time scales of electrons in the vicin- 
ity of the stars (synchrotron and IC p rocesses), the standard 
formulae (|Blumenthal fc Gouldfl97Cl ) are applied. The com- 
parison of the acceleration and the synchrotron energy loss 
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Table 2. Parameters of the nebulae around binary systems and energies of particles which produce 7-rays: the parameter 77 defining the 
shock structure, the part of the hemisphere, A m ; x , in which the mixed wind is confined, the wind mixing distance -D m i x , the magnetic 
field -B m i x at D mlx , and the bulk speed of the mixed wind w env , the acceleration efficiency of electrons \ = 10 — 3 X— 3; the- optical depth 
for IC scattering of stellar photons Tjc/ t , and the electron acceleration limit due to advection, E"^~, and synchrotron E™£f energy 
losses. 



Name 


V 


^mix 


Anix (cm) 


-Bmix (G) 




W cnv (cm/s) 


X-3 


TlC/T 


E%£< (TeV) 


EZ7 (TcV) 


B1259-63 


0.01 


0.09 


2 X 10 16 


2 x 10~ 5 - 


2 x 10- 4 


3.5 x 10 8 


0.14 


0.4 


47-150 


1.7-17 


B1957+20 


1.8 


0.82 


2.3 x 10 12 


4.4 x 10~ 4 


- 0.4 


2.7 x 10 9 


8.1 


0.4 


17-460 


0.1 - 85 


J1816+4510 


1.2 


0.95 


2 x 10 12 


4 x 10" 4 - 


0.4 


2.2 x 10 9 


5.1 


3.4 


12-410 


0.05-50 


LS 5039 


0.08 


0.32 


8 x 10 13 


8 x 10~ 3 - 


8 x 10~ 2 


1.1 x 10 9 


1.3 


110 


7-24 


7.7-77 



time scales gives us the maximum energies of electrons due 
to the synchrotron energy losses, 

« 19( X - 3 D 13 /B 2 R 11 ) 1/2 TeV. (12) 

The maximum energies due to the Inverse Compton energy 
losses (in the Thomson regime) are given by, 

£ic7t « 0A( X - 3 B 2 D ia /R 11 ) 1/2 /T! TeV. (13) 

We stress that this second limit is only important in the 
Thomson regime, i.e. for electrons with energies < IOO/T4 
GeV. For larger energies, due to the falling of the Klein- 
Nishina cross-section, the limit is relaxed. The values of the 
maximum energies of electrons due to their escape with the 
expansion of the nebula and due to the synchrotron energy 
losses are shown for the considered binary systems in Ta- 
ble. 2. We conclude that electrons can be accelerated to en- 
ergies estimated by the lower value between i?^ n x and E 1 ^?. 
Note that these maximum energies are close to those values 
expected for electrons injected from the inner pulsar mag- 
netosphere. Since it is not clear which mechanism accelerate 
electron in the vicinity of the binary system, we consider 
two limiting cases: (1) electrons are injected into the nebula 
surrounding the binary system with a power law spectrum 
extending between Ef nl = m c -yV nl and £™ ax , or (2) they 
are injected with monoenergetic energies i?™ ln . 

Relativistic hadrons are also expected to be accelerated 
in such nebulae. The maximum energies of these hadrons 
should be also limitted by their advection from the nebula. 
So then, they are given by the same formula as for electrons 
(Eg. Ill[) . We estimate the optical depth for the interaction 
of these relativistic hadrons with the matter of the mixed 
pulsar-stellar wind. It is assumed that hadrons are captured 
by the magnetic field of the mixed wind above -Dmix- They 
are advected with the velocity of the wind 

Wenv ■ 

In such a 

case the optical depth is equal to 

t pp = / n H <7 pp dD Ri3x 10- 3 ——L, (14) 

where n H = M/(47rD 2 w onv m p ) « 3 x 10 7 M- 7 /(Dl 3 wg) 
cm~ 3 is the density of the matter in the mixed pulsar- 
stellar wind as a function of the distance from the star, 
o"pp = 3 x 10 -26 cm 2 is the cross section for proton-proton 
interaction, and m p is the proton mass. For the parameters 
of the binary systems considered above (see Tables 1 and 
2), the optical depths for hadronic interactions are clearly 
lower than unity. Therefore, we conclude that relativistic 
hadrons are not able to find enough target mass within the 
considered nebulae around binary systems to efficiently lose 



energy on hadronic processes. These protons mainly escape 
from the binary system and contribute to the galactic cos- 
mic rays. Note that the stellar wind (and also mixed pulsar- 
stellar wind) can be quite inhomogeneous. However, rela- 
tivistic protons have to be very efficiently captured within 
the regions with density effectively enhanced by a factor of 
~10 5-6 , in respect to average densities estimated above for 
the considered binaries, in order to provide their efficient in- 
teraction with the matter. Such efficient capturing seems to 
be rather unlikely in the stellar winds since the filling factor 
of the stellar wind with such dense regions should be very 
low. 



5 GAMMA-RAY PRODUCTION IN NEBULAE 

We calculate the 7-ray spectra produced by relativistic elec- 
trons within the nebulae in the Inverse Compton scattering 
of the stellar radiation and the Microwave Background Ra- 
diation. We assume that relativistic electrons are isotropized 
in the mixed pulsar-stellar wind at the distances above -Dmix- 
The maximum energies of the electrons are considered in the 
range which is consistent with the estimations reported in 
the last two columns in Table 2. The electrons are captured 
by the magnetic field of the mixed winds and move grad- 
ually with the velocity of these winds. We also include the 
synchrotron energy losses of the electrons in the magnetic 
field of the winds with the values equal to B m i x at the dis- 
tance Dmix- This magnetic field drops with the distance, 
D, from the binary system as oc 1/D. For example, in the 
case of J1816+4510, we apply the value -B m ; x = 10~ 2 G. 
For this value, the energy density of stellar radiation at the 
distance D m ix overcomes the energy density of the magnetic 
field by a factor of ~200. Therefore, the synchrotron energy 
losses can only start to be important at energies above a 
few hundred GeV, where the IC emission is suppressed by 
the Klein-Nishina cross-section provided that we consider 
the upper range of the magnetic fields (equal to ~ 0.4 G, 
see Table 2). The relative importance of the synchrotron en- 
ergy losses in respect to the IC energy losses, in the case of 
LS 5039 is similar to that estimated for J1816+4510. The 
surface temperature of the companion star in LS 5039 is a 
factor of two higher but the distance of the mixing region 
from the star (in units of stellar radius) is a factor of 4 larger. 
Therefore the radiation field in the case of both sources at 
the distance D m ix are comparable. Note however, that the 
time scale for the IC scattering is a factor of ~40 larger for 
LS 5039 which affects the optical depth for electrons (see 
L'mix and Tig/t m Table 2). 
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Figure 2. 7-ray SEDs expected from the nebula surrounding the Redback millisecond binary system J1816+4510 (on the left) and 
the binary system LS 5039 supposed to contain a classical pulsar (on the right). The upper panels show the results for monoenergetic 
injection of electrons with energy: 10 2 GeV (dotted), 10 3 GeV (dashed), and 10 4 GeV (solid). The middle panels show the results for a 
power-law injection with spectral index of — 2.1 between F m i n = 10 GeV and F max — 3 x 10 2 GeV (dotted), 3 x 10 3 GeV (dashed), and 
3 X 10 4 GeV (solid). The bottom panels show the results for a power-law injection with spectral index of —3 in energy ranges: 10 — 3 X 10 2 
GeV (for J1816+4510) and 100 - 3 x 10 3 GeV (for LS 5039) (dotted), 10 - 3 x 10 3 GeV (for J1816+4510) and 10 - 3 x 10 4 GeV (for 
LS 5039) (dashed), 10 - 3 X 10 4 GeV (for J1816+4510) and 100 - 3 X 10 4 GeV (for LS 5039) (solid). It is assumed that A mix ■ e = 10% 
of the pulsar wind power is converted to relativistic electrons. In the case of J1816+4510 most of the relativistic electrons in the pulsar 
wind escapes from the binary system but in the case of LS 5039 only a part of the electrons in the w ind, estimated on A j 3 " 1 ~ 0.2, 
can escape from the binary system. For the comparison, the 100 hr differential sensitivity of MAGIC llAleksic et al.ll2012l) and CTA 
llBernlohr et alj|2012l) are also marked by the gray solid an d dotted curve respecti vely. Full and empty data circles show the SED of 
LS 5039 measured by H.E.S.S. in 2 different orbital phases d Aharo nian et alj|2006l ) and the average spectrum from 30 months of the 
observations with Ferrm-LAT is shown with triangles l|Hadaschll2012t) . The thin black line in the left panels show the spectral fit to the 
Ferrm-LAT observations of J1816+4510 l lKaplan et al.ll2012l) . 



We apply the Monte Carlo method in order to obtain 
the 7-ray spectra produced by electrons which IC scatter 
stellar photons. These electrons are advected from the bi- 
nary system with the velocity of the mixed pulsar-stellar 
wind w cnv . The turbulent magnetic field in the mixed pulsar- 
stellar wind isotropizes the directions of electrons. There- 
fore, in the reference frame of electrons the stellar radiation 
from the companion star is isotropic. We take into account 
that electrons propagate in the stellar radiation field and 
the magnetic field which decrease with the distance from 
the binary system. The 7-ray spectra are calculated apply- 
ing the formulae for the IC spectra in the general case from 
Blumenthal & Gould (1970). 



Those two binary systems, J1816+4510 and LS 5039, 
are selected for more detailed calculations. They seem to be 
the best candidates for the production of detectable 7-ray 
fluxes from the nebulae due to the luminous stellar compan- 
ions and relatively compact nebulae. We assume that rela- 
tivistic electrons take e = 10% of the rotational power lost 
by the pulsars. In the case of the monoenergetic injection 
of electrons from the pulsar magnetosphere in the binary 
system LS 5039, only electrons, moving within the solid an- 
gle directed outward from the companion star, AQjq, can 
escape from the inner region of the binary system without 
significant energy losses on the IC process. These electrons 
escape in the outward direction (in respect to the star) since 
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they move in the highly anisotropic stellar radiation field. 
Based on the previous calculations of the radiation processes 
within the binary system LS 5039 (Bednarek 2006), we esti- 
mate this part of the hemisphere on Ajq = Afi I c n /4vr ~ 0.2. 
Thus, only 20% of the relativistic electrons escape without 
significant energy losses to the outer regions of the binary 
system LS 5039. This additional factor has been included 
in our calculations for LS 5039. In the case of J1816+4510, 
electrons moving rectilinear through the binary system do 
not lose efficiently energy on the IC process due to lower 
surface temperature and smaller radius of the companion 
star. Therefore, we assume A^ n ~ 1 for J1816+4510. 

The 7-ray spectra, expected from the nebulae around 
these two binary systems, are shown in Fig. 2 for the case 
of the monoenergetic spectrum of electrons (upper figures) 
and the power law spectra of electrons which were ejected 
with the spectral indices equal to 2.1 (middle figures) and 
3 (bottom figures). The parameters defining these spectra 
are reported in the figure caption. The maximum energies of 
electrons are consistent with the limits derived for the case of 
dominant synchrotron energy losses or advection time scale 
with the mixed winds (Table 2). 

The 7-ray spectra, expected at the observer, are com- 
pared with the sensitivity of the present generation (e.g. 
MAGIC) and future (Cherenkov Telescope Array, CTA) 
imaging atmospheric Cherenkov telescopes (IACTs) . Detec- 
tion of the considered here TeV 7-ray emission from the 
Redback millisecond binary JI8I6+4510 with the present 
instruments is unlikely. However, with a deep exposure of 
100 hr, such emission could be detected by CTA as long as 
the energy conversion efficiency, e, is at least 10%. The 7- 
ray spectrum obtained from the Fermi-LAT observations of 
that source does not constrain the product of part of the 
hemisphere subtracted by the mixed wind and the energy 
conversion efficiency for this source, i.e. A m i x • e. 

The 7-ray spectra expected from the LS 5039 are ~2 
orders of magnitude above the expected sensitivity of CTA. 
With the assumed energy conversion efficiency from the pul- 
sar to relativistic electrons (10%), such emission should be 
clearly detected even by the presently operating Cherenkov 
telescopes. The TeV emission observed from LS 5039 by the 
H.E.S.S. telescopes is modulated with the orbital period. On 
the other hand, the emission expected from the presented 
here model should contribute as a constant component in 
the whole emission of the system. Therefore, we can use the 
low state observed by H.E.S.S. telescopes to constrain the 
parameters which determine the energy transferred to the 
TeV electrons which escape from the binary system. This 
power depends on the product of a part of the hemisphere 
in which the mixed stellar-pulsar wind is confined, A m i x and 
the energy conversion efficiency, e. Derived upper limits on 
A m i x ■ e are reported in Table. 3 for different spectra of 
electrons in the case of LS 5039. For the broad range of 
the model parameters, A m ; x • e is severely constrained to be 
« 1%. 

It is believed that the energy conversion efficiency from 
the pulsar to relativistic electrons, in the inner pulsar mag- 
netosphere and in the shock acceleration scenario, is of the 
order of 10%. Then, the part of the hemisphere in which the 
mixed wind is confined should be also below A m i x ~ 10%. 
As shown in Sect. 2, this factor depends on the half opening 
angle of the shock within the binary system which is de- 



termined by the value of the parameter rj. For the value of 
rj, calculated for the supposed shock structure in LS 5039, 
A m i x is larger than required above (see Table 2). In fact, 
different phenomena can effect the estimated value of A m i x . 
At first, the value of rj might be much larger than ~0.08, 
i.e the stellar wind may not effectively confine the pulsar 
wind. In fact the winds around Be type massive stars are 
expected to be aspherical with the dense and slow equato- 
rial wind and the fast and rare poloidal wind. If the pulsar 
is mainly immersed in the poloidal wind, than the shock 
st ructure might even bend around the sta r (see calculations 
in ISierpowska-Bartosik fc Bednarek! 120081 ). i.e. the value of 
rj is larger than unity. In such case only a part of the pulsar 
wind can mix with the stellar wind. Second, the electrons 
could be injected anisotropically from the pulsar or acceler- 
ated anisotropically in the pulsar wind. In fact, the pulsar 
winds are expected to be highly anisotropic with the dom- 
inant equatorial compo nent fe.g. iBogovalov fc Khangoulianl 
l2002l ; IVolpi et alj|2008h . Such aspherical pulsar wind will in- 
troduce significant complications to the discussed scenario. 
More energy in relativistic electrons will be directed to the 
equatorial regions of the binary system resulting in much 
larger effective value of rj (even larger than unity). However 
the solid angle, in which mixed stellar-pulsar wind is con- 
fined, will become lower. Simple estimations, which base on 
the formulae in Sect. 2, show that the increase of significant 
effective power of the wind in the equatorial direction results 
in more power transferred in the form of relativistic electrons 
to the nebula. However, in such a case a significant part of 
the particles, accelerated by the pulsar towards the compan- 
ion star, could not be able to escape from the binary system 
due to efficient IC interactions with the stellar radiation. 
Therefore, the final effect of the anisotropic pulsar wind on 
the energy transferred to nebula in the form of relativistic 
particles is difficult to estimate without more detailed calcu- 
lations. The third possibility is that more electrons could be 
accelerated in the general direction towards the companion 
star. This might happen in the case of a significant collima- 
tion of the relativistic particles in the pulsar winds by the 
shock structure (e.g. see lDubus et al. I l2010l : IBogovalov et al.l 
120081 . Then, electrons lose efficiently energy on IC process 
before they manage to escape into the mixed wind region. 

We conclude that the derived limits on the parameter 
A m ; x • e, determining the high energy emission from the neb- 
ula around LS 5039, indicate on either different parameters 
of the stellar wind than usually expected in the literature, or 
on the importance of anisotropic stellar and pulsar winds, or 
on the anisotropic injection/acceleration of relativistic elec- 
trons within the binary system. 



6 CONCLUSIONS 

We considered nebulae around binary systems which con- 
tain rotation powered pulsars and hot companion stars with 
a relatively large mass loss rate. Such nebulae differ signifi- 
cantly from the pulsar wind nebulae around isolated rotation 
powered pulsars since a part of (or whole) the pulsar wind 
becomes loaded with the mass from the companion star. As 
a result, the pulsar wind relativistic plasma slows down and 
mixes with the matter of the stellar wind. The mixed pulsar- 
stellar wind moves from the binary system in the outward 
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Table 3. The upper limits on the product of the part of the hemisphere overtaken by the mixed pulsar-stellar wind (A m ; x ) and the 
energy conversion efficiency (e) from the pulsar to relativistic electrons for the binary system LS 5039 obtained from comparison with 
H.E.S.S. and Fermi-LAT spectral measurements 







Monoenergetic 




Power-law, a = —2.6 


Power-law, a = —3 




Smax 
-^min 


[GeV] 
[GeV] 


10 2 10 :i 
10 2 10 3 


10 4 
10 4 


10 2.5 10 3.5 10 4.5 
10 10 10 


10 3.5 1Q 4.5 

100 10 


10 4.5 

100 


^mix 


■ € 


< 0.3% < 0.05% 


< 0.03% 


< 0.3% < 0.3% < 0.2% 


< 0.4% < 0.3% 


< 0.4% 



direction with a reduced velocity determined by the power 
of the pulsar wind and the mass of the stellar wind. We ar- 
gue that the distance at which these two winds have to mix 
completely is determined by the orbital period of the binary 
system. This mixing process can concern the whole pulsar 
wind (for r\ < 1 the stellar wind pressure dominates over the 
pulsar wind pressure) or only a part of the pulsar wind for 
77 > 1. In the latter case, the rest of the pulsar wind (along 
polar regions of the binary system) expands freely interact- 
ing with the interstellar medium at large distances from the 
binary system as observed in nebulae around isolated pul- 
sars. 

Due to the efficient mixing, relativistic electrons are 
isotropized in the reference frame of the mixed wind rela- 
tively close to the binary system. Therefore, they can effi- 
ciently comptonize thermal radiation from the companion 
star. This is in contrary to the nebulae around isolated pul- 
sars where the inner part of the nebula is non-radiative due 
to the almost radial motion of the plasma from the pulsar. 

We calculate the 7-ray spectra produced by leptons in 
the IC scattering process in the case of the mono-energetic 
particles (leaving the light cylinder radius of the pulsar) and 
in the case of the power law spectra of leptons formed as 
a result of re-acceleration at the shocks produced in the 
pulsar-stellar wind collisions. Note that this 7-ray emission 
should be steady, i.e. independent on the phase of the binary 
system. As an example, we show the 7-ray spectra for the 
supposed pulsar in the binary system LS 5039 and for the 
recently discovered Redback type millisecond binary system 
J1816+4510, containing relatively hot companion star. It is 
found that if electrons take a part of the pulsar wind energy 
loss rate defined by the product of the part of the hemi- 
sphere, in which the mixed pulsar-stellar wind expands, and 
for the acceleration efficiency of electrons AQ,-e equal to 10%, 
than the TeV 7-ray emission from the millisecond binary 
system is comparable to the sensitivity of the Cherenkov 
Telescope Array (CTA). However, in the case of LS 5039, 
electrons with such value of Afl ■ e should produce 7-ray 
emission on the level clearly above the sensitivity of both 
the present Cherenkov telescopes and CTA. By comparing 
the 7-ray spectra predicted for LS 5039 with the TeV 7-ray 
emission in the low state established by the H.E.S.S. tele- 
scopes, and with the average Fermi-LAT observations, we 
put strong limits on A m i x • e < 1%. Such low limit suggest 
that either the parameters of the stellar wind are less ex- 
treme than considered up to now in the literature or the 
injection of electrons from the inner pulsar magnetosphere 
and/or accelerated at the shock is highly anisotropic. Most 
of relativistic electrons might be directed towards the com- 



panion star and lose efficiently energy already within the 
binary system. 

It seems that the TeV 7-ray binary containing classi- 
cal radio pulsar PSR B1259-63 should be also a good can- 
didate source for the application of the above considered 
model. However, this pulsar is on a very elongated orbit 
around the companion star. Therefore, the nebula around 
this binary system is expected to have complicated, non- 
spherically symmetric shape. In this paper we provided cal- 
culations only for the binary systems which are expected to 
produce nebulae which can be more or less approximated 
by a symmetric shape. The more complicated case of the 
nebula around binary system, such as that containing PSR 
B1259-63, will be discussed in detail in a future paper. 
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